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Brittle materials exhibit sharp dynamical fractures when meeting Griffith's criterion, whereas 
ductile materials blunt a sharp crack by plastic responses. Upon continuous pulling ductile materials 
exhibit a necking instability which is dominated by a plastic flow. Usually one discusses the brittle 
to ductile transition as a function of increasing temperature. We introduce an athermal brittle to 
ductile transition as a function of the cut-off length of the inter-particle potential. On the basis 
of extensive numerical simulations of the response to pulling the material boundaries at a constant 
speed we offer an explanation of the onset of ductility via the increase in the density of plastic modes 
as a function of the potential cutoff length. Finally we can resolve an old riddle: in experiments 
brittle materials can be strained under grip boundary conditions, and exhibit a dynamic crack when 
cut with a sufficiently long initial slot. Mysteriously, in molecular dynamics simulations it appeared 
that cracks refused to propagate dynamically under grip boundary conditions, and continuous pulling 
was necessary to achieve fracture. We argue that this mystery is removed when one understands 
the distinction between brittle and ductile athermal amorphous materials. 



I. INTRODUCTION 

Consider the following numerical experiment, which is 
also exhibited in Fig. [1] and in the movies available in 
[l| . A piece of amorphous solid at rest is contained in the 
square (pink) box as plotted. The material is prepared 
by quenching a binary liquid of point particles (see po- 
tentials below), and removing the periodic boundaries. 
The resulting system has side walls consisting of pinned 
particles (identical to the bulk particles, blue and brown 
in Fig. [T]) forming a rigid body which is pulled, the upper 
and lower boundaries are free. The walls are moved by 
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FIG. 1: Color online: The result of pulling quasi-statically 
the side boundaries (the particles in the side walls are held 
frozen). Upper panel: the material is "brittle" and the pulling 
results in a clean crack. Lower panel: the material is "duc- 
tile" and the pulling induces a plastic flow with a necking 
instability. Note the large difference in the relative increase 
in system length compared to the original system size in the 
two panels. The number of particles TV = 800 



an increment equal in size and opposite in sign. 

In the upper panel the material is brittle, and cannot 
deform much before it fractures. In the lower panel the 
material is ductile, and the distance between the walls in- 
creases by about a factor of two before the plastic defor- 
mation that is induced creates a thin neck that is eventu- 
ally broken. The huge difference in response is obtained 
with materials in the two panels that are identical in ev- 
ery respect except for the cutoff length of the attractive 
part of their inter-particle potentials, see Fig. [2] The aim 
of this paper is to explore the intrinsic relation between 
the microscopic properties of the inter-particle potential 
and the brittle to ductile transition. We explicitly do not 
consider in this paper the interesting effects of tempera- 
ture, aspect ratio, system size and dimensionality, which 
of course are also important in determining a material's 
brittleness. 

One of our main motivations in this paper is a mystery 
that had clouded the fracture literature for a number of 
years now. The issue is crack propagation in amorphous 
solids under grip boundary conditions. In experiments 
one can pull the boundaries of a given sample and grip 
them, and then cut a slot at the edge of the sample [2j. 
In brittle materials this leads to a fast development of 
the a dynamic fracture (when the length of the initial slot 
exceeds the Griffith's length jU ) . For some unknown rea- 
son one always failed to repeat such experiments in direct 
molecular dynamic simulations, leading recently to some 
claims that only by defining 'bonds' between particles 
and letting these bonds 'break' under some criteria one 
can sustain a dynamical brittle fracture Obviously 
this is not a satisfactory answer; here we argue that the 
mystery is removed when one understands what makes a 
material brittle or ductile, and in truly brittle materials 
there is no problem in seeing a dynamical fracture under 
grip boundary conditions. 

In this paper we want to understand how the response 
of an amorphous system subjected to uniaxial load de- 
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FIG. 2: Color online: Upper panel: Examples of three dif- 
ferent potentials with r co = 1.2, 1.8 and 2.4 used in our nu- 
merical experiments. The first and the third correspond to 
the experiment described in Fig. [1] Lower panel: the pair 
radial distribution function g(r) for the three corresponding 
potentials in the upper panel. Note that the second shell of 
neighbors begins around r = 1.5 



The structure of this paper is as follows: In Sect. |TT] 
we describe the glassy system used and how we control 
the range of the inter-particle potential. The numerical 
experiments and their results are described in Sect. IIIII 
In Sect. IIVI we relate and explain the findings of Sect. 
IIIII to the density of states of the glassy materials, ex- 
plaining that with increasing the range of the potential 
one changes the weight of the contribution of the plastic 
modes in the density of states. This correlates very well 
with the tendency towards ductility as seen in the numer- 
ical experiments. Finally in Sect|V]we focus on the riddle 
of brittle fractures and show that there is no riddle at all 
- one simply needs to make the material brittle enough 
to allow brittle fractures to run spontaneously under grip 
boundary conditions. In Sect. I VII we conclude the paper 
and reiterate the dramatic effect of the inter-particle's 
potential properties on the behavior of a material under 
load 0. 

II. SYSTEM AND METHODS 



pends on the parameters of the potential, and in partic- 
ular the cut-off length r co . This will provide us with a 
fair insight into the athermal brittle to ductile transition. 
We will relate the degree of ductility to the growing im- 
portance of the contribution of "plastic modes" to the 
density of states of the material. 



For the numerical experiments we employ a generic 
glass former in 2-dimensions in the form of a binary mix- 
ture whose amount of bi-dispersity of 'small' and 'large' 
particles was chosen to avoid any crystallization. In fact 
the particles interact by inter-particle potentials as shown 
in Fig. [2j with the analytic form 
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Here r m - m /Xij is the length where the potential attain 
it's minimum, and r co /Xij is the cut-off length for which 
the potential vanishes. The coefficients a, b and c^i are 
chosen such that the repulsive and attractive parts of 
the potential are continuous with two derivatives at the 
potential minimum and the potential goes to zero contin- 
uously at r co /Xij with two continuous derivatives as well. 
The interaction length-scale Xij between any two parti- 
cles i and j is Xij — 1.0A, Xij — 1.18A and Ay = 1.4A for 
two 'small' particles, one 'large' and one 'small' particle 
and two 'large' particle respectively. The unit of length 
A is set to be the interaction length scale of two small 
particles, e is the unit of energy and kg = 1. 

The systems employed in the pulling experiments con- 
sist of configurations which were quenched from systems 
initially equilibrated at high temperature and pressure 
(T = 1.0, P = 1.0) with periodic boundary conditions. 



The number of particles N varied between N = 800 
and TV = 67000. We used a modified Berendsen ther- 
mostat which couples a constant number of particles to 
the bath, regardless of the system size @. The systems 
are cooled to temperature T = 10~ 2 while keeping the 
pressure high, in order to avoid the creation of any holes 
in the material. Then the pressure is reduced to zero, 
P = 0.0, such that the periodic boundary conditions 
could be removed; particles forming the right and left 
walls were frozen, but the upper and lower boundaries 
were rendered free. 

From this point the system was treated differently for 
the brittle crack and the necking experiments. The brit- 
tle crack experiment starts by loading the system uniax- 
ially (with a constant velocity such that D wa ii c s /10) 
until a desired stress is reached, and then the side walls 
are held fixed. A cut is then implemented by the cancela- 



3 




FIG. 3: Color online: The dependence of r max on r co for 
the geometry of Fig. [T] (square initial conditions). Note that 
Fmax tends to saturate at a given value r^ ax which becomes 
independent on r co ; in athermal conditions it depends on the 
geometry of the sample only. 



tion of forces crossing an imaginary line of desired length 
which starts at the lower boundary. The evolution of the 
crack is simulated by molecular dynamics at T = 0.0 and 
requires no further loading of the system. 

The necking experiments were conducted in athermal 
quasi-static conditions. For those experiments the sys- 
tems were further minimized with respect to the acting 
forces on each particle, making them athermal. The right 
and left walls were then moved with a strain increment 
81 = 5 • 10~ 3 . The response we investigate is the rela- 
tive increase in the system's length T = AL/L along the 
stretching direction. The maximal value that the system 
is able to accommodate before it breaks into two parts is 
denoted as r max Q. 



III. NUMERICAL EXPERIMENTS AND 
RESULTS 

In this paper we describe only 2-dimensional simula- 
tions with a square initial geometry. In fact we have 
set up samples of height H and length L = aH, with 
a G [1,5], but we did not find any exciting aspect-ratio 
dependence beyond the obvious. For the square initial ge- 
ometry the athermal response was characterized by F max . 
In Fig. [3] we exhibit the dependence of r max on r co for 
square systems (a = 1). It is sufficient to limit the range 
of the potential cut-off to 
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The lower limit for r co is dictated by the position of the 
first peak in g(r) cf. Fig. flower panel. For r co < 1.2 
the potential gets shorter than the first shell of neigh- 
bors. The upper limit is determined by the convergence 
of our 1-parameter family of potentials to the asymptotic 
Lennard- Jones potential. Within this range of potentials 
(r max ) was evaluated by averaging over 4500 samples. It 
is found to increase smoothly with r co reaching a limiting 



value r^ ax . When the system reaches the limiting value 
r™ ax we refer to it as "maximally ductile" . 

The qualitative understanding of the result shown in 
Fig. [3] is quite straightforward. In the limit r co — > r m ; n 
any cavity that forms whose dimensions is larger than r co 
will not heal since there are no attractive forces operating 
across it. Rather, it will propagate rapidly to span the 
system in a form of a brittle fracture. With r co increas- 
ing, a small cavity will still be held by attractive forces 
across it. The higher is r co the larger is the number of 
neighbors that each particles feels attracted to, and one 
needs to create a very large cavity indeed to overcome 
these attractive forces. On the other hand when r co in- 
creases beyond the second shell of the binary correlation 
function, not much can change. With our choice of po- 
tentials, in the limit r co — > oo we simply converge to 
the infinite-range Lennard- Jones potential which has a 
negligible attractive force beyond, say, r/X = 3. This 
qualitative understanding is further supported by exam- 
ining the contribution of the plastic modes to the density 
of states, as we do next. 



IV. THE CONTRIBUTION OF PLASTIC 
MODES 

It is well known that in amorphous solids the density 
of states deviates from the classical Debye formula which 
was developed for a perfect elastic solid. The deviations 
were referred to as the 'Boson peak' but very little reliable 
analytic statements were made about the form of this 
deviation Q. In recent work of the present group it was 
proposed that one can provide a useful analytic form of 
the density of states (the HKLP model @) which will be 
particularly useful at the low end tail of frequencies. In 
terms of the eigenvalues of the Hessian matrix Aj = uif 
where is the i'th normal mode frequency the proposed 
model is 
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where the first term on the RHS is the Debye formula 
in which N is the number of particles, d the space di- 
mension and Ad is the Debye cutoff frequency which in 
2-dimensions reads Xr> = 87r/i with /i being the shear 
modulus. The second term is our proposed addition to 
the model density of states that should be appropriate for 
all generic amorphous solids in the limit A -C Ajj. This 
term contains contributions from eigenvalues Xi that de- 
pend on the external mechanical loading, and go to zero 
with increasing the load, harking the appearance of plas- 
tic instabilities 9]. The exponent 9 in this plastic con- 
tribution is not universal, it depends on the amount of 
disorder in the glassy system, and it was shown to be an 
important characteristic of the mechanical properties of 
amorphous solids [9j. 

In numerical calculations it is easier to evaluate the 
cumulative distribution C(A) = 



/ A P(X)dX, whose form 
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FIG. 4: Color online: the cumulative density of states from 
which the Debye contribution was subtracted. The remaining 
plastic contribution is plotted in log-log coordinates such that 
the intercept is B/(l + 8) and the slope is 8. 




FIG. 5: Color online: the value of the weight B in Eq. [4] as a 
function of r co . Inset: r max vs. B 2 . 
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In the present context we wish to understand the re- 
lation between the tendency towards ductility as r co in- 
creases to the amount of plasticity as expressed, say, by 
the value of the weight coefficient B and the exponent 9. 
To this aim we have computed the cumulative density of 

states C \ an d subtracted from it the exactly known 

Debye contribution. The resulting plots are shown in 
Fig. [4] We see that the model[3]is faithfully representing 
the data for small values of r co , whereas for the largest 
ones the amount of plastic modes increases beyond the 
first term contribution in the model. This is however not 
causing any difficulty in estimating the weight B (the 
intercept in Fig. @| which is shown in Fig. [5] 

It is quite obvious that our T max dependence on 7- co 
(cf. Fig. [3]) is qualitatively similar to the dependence of 
the weight B on r co . It is therefore natural to plot r max 
as a function of B. A convincing linear dependence can 
be obtained by plotting r max as a function of B 2 as is 
shown in Fig. \5\ 
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FIG. 6: Color online: the development of a brittle crack under 
grip boundary conditions. The finite size effect appears as 
bending of the upper boundary, releasing the stress before 
the fracture hits the boundary. 



V. THE BRITTLE CRACK 

We are now able to construct a system where all the 
microscopic parameters are tuned to minimize the occur- 
rence of plastic events, a system which can be denoted 
'maximally brittle'. We expect such a system to fail un- 
der grip boundary conditions upon the application of a 
cut in one of the free boundaries. For a large value of 
r co , a small cavity can still be held by attractive forces 
across it. The longer is r co the larger is the number of 
neighbors that each particles feels attracted to, and one 
needs to create a very large cavity indeed to overcome 



these attractive forces. On the other hand, in the limit 
of r co — > r m i n any cavity that forms whose dimensions is 
larger than r co cannot heal since there are no attractive 
forces operating across it. Rather, it would propagate 
rapidly to span the system in a form of a brittle frac- 
ture. The system we used to demonstrate this concept is 
composed of N — 67000 particles with a cut-off length of 
r co = 1.20 and e = 10.0; this choice was made to increase 
the stability of the system and prevents any unwanted 
events during the loading part of the simulation. The 
experiment is run by means of molecular dynamics at 
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FIG. 7: Color online: the development of a brittle crack ap- 
pears many times to consist of the creation of a cavity ahead 
of the crack, followed by a necking process in which the cav- 
ity and the body of the crack coalesce. Shown are two such 
events, where the left panels exhibit the cavities ahead of the 
crack, and the right panels exhibit the crack after the necking 
events coalesced the crack with the cavities. 

T = 0.0. The size of the cut is 100A. The evolution of 
the crack can be appreciated in Fig. [51 

When the fracture originates from a cut, the system 
does not need to create a cavity in order to propagate a 
crack. The stress field around the tip of the cut is greatly 
enhanced by the strong curvature of the boundary, and 
once the Griffith's criterion is fulfilled the crack can run. 
Nevertheless during the propagation of the crack we ob- 
serve many times the formation of a cavity ahead of the 
crack tip followed by necking to coalesce the cavity with 
the body of the crack. The observed dynamics is pre- 
sented in Fig. which employs snapshots from the brit- 
tle crack development that is shown in Fig. |51 We note 
in passing that this mechanism of propagation of brittle 
cracks was seen experimentally in Ref. [lfj. In Ref. 
this same mechanism was proposed as the reason for the 
roughening of the fracture surface. This roughening is 



seen to the bare eye in Fig. [6j its quantitative assess- 
ment and its relation to the mechanism in question will 
be considered in a future publication. 



VI. SUMMARY AND CONCLUSIONS 

The main result of this paper is the dramatic change in 
material mechanical properties that results from chang- 
ing the inter-particle potential interaction length. We in- 
troduced the athermal brittle to ductile transition, show- 
ing that even at T = we can go all the way from brittle 
to ductile response only by changing the cutoff length r co . 
There is an interesting analogy between increasing the 
cutoff length and raising the temperature; in both case 
the system becomes 'softer'. This softening can be quan- 
titatively described by our model of density of states of 
amorphous solids, in which the coefficient B measures the 
increasing weight of plastic modes, which are those eigen- 
functions of the Hessian matrix whose eigenvalues vanish 
in saddle-node bifurcations when the system is strained 
@. Indeed, we found a linear dependence between the 
degree of lengthening of the system before breaking and 
B 2 . Finally, the new understanding of brittleness vs. 
ductility allowed us to disperse a long standing conun- 
drum - how to get a brittle crack running under grip 
boundary conditions. It turns out that previous attempts 
failed simply because the system chosen were not brittle 
enough. Building a 'maximally brittle' material results 
in a nice running crack under grip boundary conditions 
as seen in many experiments. 
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